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PACS 05.50.+q - Lattice theory and statistics (Ising, Potts, etc.) 
PACS 05 . 70 . Fh - Phase transitions: general studies 
PACS 64 . 60 . ah - Percolation 

Abstract. - We present a study of phase transitions of the mean-field Potts model at (inverse) 
temperature (3, in presence of an external field h. Both thermodynamic and topological aspects of 
these transitions are considered. For the first aspect we complement previous results and give an 
explicit equation of the thermodynamic transition line in the (3—h plane as well as the magnitude 
of the jump of the magnetization (for q ^ 3). The signature of the latter aspect is characterized 
here by the presence or not of a giant component in the clusters of a Fortuin-Kasteleyn type 
representation of the model. We give the equation of the Kertesz line separating (in the /3-h 
plane) the two behaviours. As a result, we get that this line exhibits, as soon as q 3, a very 
interesting cusp where it separates from the thermodynamic transition line. 
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Introduction . — In [1], Kertesz pointed out that 
a very interesting phenomenon arises in the Ising model 
subject to an external field. The so-called Coniglio-Klcin 
droplets [2] associated to Fortuin-Kasteleyn clusters [3] 
have a whole percolation transition line extending from 
the Curie point to infinite fields. This seems, at first sight 
in contradiction with the fact that on the other hand, ther- 
modynamics quantities do not have any singularities in 
any of their derivatives (with respect to the temperature 
or the field) as soon as the field is non-zero. But, as 
Kertesz already mentioned: 

However, we emphasize that the suggested pic- 
ture is not in contradiction with the non ex- 
istence of singularities in the thermodynamic 
quantities because the total free energy remains 
analytic. 

Indeed in such a model, the criticality can be specified in 
two different ways: the thermodynamic criticality associ- 
ated to the thermodynamic limit of the bulk free energy 
and the geometric criticality associated in the same limit 
to another lower order free energy. 

Since then, the Kertesz line has remained the subject 
of interests over the years for various models of statistical 



mechanics. In particular, it has been considered recently 
in [4] within the Potts' model on the regular lattice Z d . 
There it was found, that, in dimension d = 2, the whole 
thermodynamic first order transition line (when such first 
order behavior is present in the system) coincides with 
the Kertesz line. The latter is also first order in the cor- 
responding range of values of temperature and field. 

The aim of this letter is to present a study of such a 
model when the underlying lattice is the complete graph 
with n vertices. In this case the model is called mean-field 
or Curie- Weiss Potts model. 

The model. — To introduce the Potts model on the 
complete graph we attach to the sites i = l,...,n of 
the graph, spin variables <7j that take values in the set 
{1, . . . ,<?}. The model at temperature T = 1//3 and sub- 
ject to an external field H = h/(3 is then defined by the 
Gibbs measure 
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over spins configurations cr. Here, Zp tts denotes the par- 
tition function, the indices i, j runs over the set {1, . . . , n}, 
and S is the Kroncckcr symbol. The critical (thermody- 
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namic) behaviour of this model is well known [5], [6], [7], 
and mainly governed by the mean field equation 



-/3 s + In 



1 + (q - l)s 



1 - s 



Namely, when h = 0, there exists a threshold value 



2 for 
.0-1 



q-2 



\n(q — 1) for q ^ 3 



such that the system exhibits at /3c a continuous transi- 
tion when q = 2, and a first order transition when q 3. 
When h > 0, no transition occurs if q = 2, while as soon 
as q ^ 3 a first order transition line appears for which the 
microcanonical free energy of the model has two minima 
associated with two different solutions of the mean field 
equation [5], [7]. 

To study the behaviour of clusters previously men- 
tioned, we turn to the Edwards-Sokal joint measure [8] 
given by 



(2) 
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where the edges variables r)ij belong to {0, 1}. 

Notice that when h = oo, all 04 = 1. This means that we 
open edges (rjij = 1) with probability (w.p.) p = 1 — e _/3 /™ 
and close edges (rjij = 0) w.p. e"' 3 /" . This is nothing else 
than the well known Erdos-Renyi random graph Q(n,p) 
[9] . This random graph is known to exhibit a (topological) 
transition at [3 = 1 such that with probability tending to 
1 as n — > oo: 

i. for (3 < 1 all components of open edges are at most of 
order Inn. 

ii. at P = 1 a giant component of order n 2 / 3 appears. 

iii. for (3 > 1 this giant component becomes of order s*n 
where s* is the largest root of the mean field equation 
(2) with q = 1 and h = 0. 

We refer the reader to [10], [11] for detailed discussions and 
proofs (see also [12] for a new approach of this transition). 

Notice also that, on the other hand, when h — 0, the 
marginal of the ES measure over the edges variables is the 
random cluster model: 

mc(v) = J- n <r(/Vn)(1 ~" w) ( 1 - e-^q CM (6) 



is given in [13] (see also [14], [15]). There, it is shown that 
with the threshold value (3 c q \ given by (3) for < q ^ 2 
and by (4) for q > 2, then with probability tending to 1 
as n — > oo: 

a) if j3 < [3c q \ the largest component (of open edges) of 
/xrc is of order Inn. 

b) if f3 > /3c q \ the largest component of /irc is of order 
son, where so > is the largest root of the mean field 
equation (2) with h = 0. 

c) if (3 = (3 ( c q) and < q < 2, urc has largest component 
of order n 2 / 3 . 

d) if (3 = [3^ and q > 2, /xrc is either as in a) or as in b). 

According to these results, it is natural to expect for 
model hes (and its marginal over the edges variables), 
a Kertesz line hxiP), interpolating betwen /jk-(I) = oo 
and hfc{Pc) = 0, that signs the emergence of a giant 
component in the open edges of the configuration ry. 

To study this line, it is convenient to consider a col- 
ored version of the Edwards-Sokal formulation. Namely, 
whenever the endpoints ij of a given edge are occupied 
by a spin of color a, we label the variable rjij with a su- 
perscript indicating the color. For a spin configuration 
er, let n a (er) be the number of sites occupied by spins of 
color a = 1, . . . , q. We then relabel the variables by 
r/jy : k < I S {1, . . . ,n (er)}. This means that, for any 
pairs i < j G {1, . . . , n a (cr)}, we open an edge r)£j = 1 

w.p. 1 — e~^/ n , and close this edge 77^- = w.p. e~Pl n . 
All the other edges between two sites of different colors 
are closed. The resulting measure becomes 
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i<je{l,...,n a (o-)} 

By summing over the spins variables, we get the following 
Colorcd-Random-Clustcr model 
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where x a = n a /n are the densities of the colors and 



Z(ni, ...,n q ) 
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where C(rj) denotes the number of connected components 
(including isolated sites) of open edges of the configuration 
77. For q = 1, this model again reduces to Q(n,p) withp as This colored random cluster model can be thought as fol- 
before. A refined study of the random cluster model (6) lows. Given a partition m, . . . , n q of n, we have q classical 
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Erdos-Renyi random graphs where the edges are open w.p. 
1 — e~ /3Xa / r ' a . The asymptotic behavior of the partition as 
n — > oo will be determined by minimizing the free energy 

f(xi,...,x q ) = - lim -\nZ{n\,...,n q ) (10) 

n— *oo n 

where the limit is taken in such a way that n a /n — » x a . 

Notice that, up to a normalizing factor, Z(n\, . . . , n q ) is 
the microcanonical partition function of the Potts model 
(1) restricted to configurations such that the number of 
spins of color a is fixed to n a . This will give the thermo- 
dynamic behavior of the system and will determine the 
densities x a . Once this is done the topological properties 
can be analysed by using the known properties of classi- 
cal random graphs. Due to the presence of the field or 
using the symmetry of colors for vanishing field, we will 
pick up the first color. Then, the threshold value of the 
(topological) transition will be given by 



1 



(11) 



For A = /3xi > 1 a giant component will appears, with 
size of order & Xl n = 0px 1 xin where 9\ is the largest root 
of the mean field equation 



xe 



ln(l - 6) 
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Results. — Let us first consider the case q = 2. 

When h = 0, the mean field equation (2) has a unique 
solution so = for /3 ^ 2, and a unique solution sq > 
for (3 > 2. Inserting sq into (17), we obtain that the 



giant component appears at (3 K = 2, getting thus 0%- 

P^ 2 \ There, this component is of order n 2 / 3 . For /3 > 2, 
Pxi(sq) > 1, and we are in the supercritical regime with a 
giant component of order Q xl ( so )n. For (3 < 2, (3x\(sq) < 
1, and we are in the subcritical regime with components 
of order at most Inn. 

When h > 0, the mean field equation has a unique so- 
lution so > for any f3 (no thermodynamic transition oc- 

(2) (2) 

curs). The transition line (3 K (h), or equivalently h y K '((3), 
is obtained by eliminating s between (2) and (17). This 
Kertesz line is given by 



,(2) 



13-2- ln(/3 - 1) 



(18) 



On this line, the largest component is of order n 2 / 3 . For 
h < hft, the largest component is of order Inn, while for 
h > h^i the largest component is of order Xl ( So \n. 

Let us then consider the case q ^ 3. We have now a 
thermodynamic transition line 



Let us briefly recall the minimization procedure for h,T\p) — —ft 
(10). By Stirling's formula, we have n\jn\ 



g-2 
2(9-1) 



+ ln(g-l), (19) 
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at the leading order in n, giving 



) = zZ Xa lnXa ~ ? zZ x * ~ hx 



(13) 



with endpoints 



(Po = 4 q —,h Q = -2 q — 



ln(q-l)), ((3^,0) (20) 



The minima of this function can be parametrised by a 
real number s, s ^ 1, such that the minimizing vector 
(x\, . . . , x q ) have the components 



Xl 



l + (q- l)s 



1 -S 
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(with arbitrary order for h = 0). 

The free energy (13) then takes the form 



1 +(<?-!> 
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9-1 



ln(l + (q- l)s) 
s )]n{l-s)-j3^^-s 2 
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The minimizers have to be found among the solutions of 
the mean field equation (2) obtained by differentiating / 
with respect to s. Note that by taking into account (14) 
the condition (11) becomes 



K - 



1 + (q ~ l)a 



1 



(17) 



As previously mcntionned, the thermodynamic transition 
has been already established in [7], for a curve with end- 
points (20). To show that the thermodynamic curve is the 
straight line (19) we observe, on hx the free energy (16) 
has the symmetry /(s) — f!3z— 
Fig. 1. 



/(fEf-a) with se[0,fEf], see 




q-2 



q-2 



2(q-l) q-1 

Fig. 1: The free energy on /it- 
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On hr the free energy has two minima s_ and s+ for 
which /(s_) = /(s_|_). They satisfy the mean field equa- 
tion (2) and thus are given by the parametric equations 



(fi,h) = (fi s± ,- 



2(9-1) 



fi s± + 111(9 - 1)) (21) 



where 



2(9 -i; 



(q - 1)(1 - s) 



and s_ G [0, 2^)J> s + e L2( g -i)> (^=TJJ 
For /3o < fi ^ /3c (or equivalently ^ h < ho) these two 
minima are distinct: s_ < s + . At (fio,ho), s_ = s + = 
and at (fi { c q) ,0), s_ = and s+ = 

Outside of the segment (19), there is only one minimum 
so- This minimum is an analytic function of fi and /i, and 
so — ► s_ as h 1 Ht, so s+ as /i J. /it- In addition, on 
/iTj as /3 increases from /?o to filf 1 (or equivalently as /i 
decreases from ho to 0), fixi(s-) strictly decreases from 
fio/2 > 1 to fie /q < 1, while fix\(s + ) strictly increases 
from fio/2 to fi c q) (q - l)/q. See below. 

We now turn to the topological behavior. 
Let (ficp, h cp ) the point of /it for which fix\{s-) = 1. This 
point is distinct from the two endpoints (20). It is given by 
the solution of the equation (/3 cp — 2)g = 2(q— 1) ln(/3 cp — 1) 
which can be solved by using the Lambert W function. 

As a consequence of the above remarks, one gets the 
following topological behaviour for the system. On the 
thermodynamic line hx- 

1. at (ficp,h C p), the largest component of /icrc is either 
of order n 2//3 or of order @ Xl ( s+ )n. 



q—2 q-2 



2. if h, 

cither of order Q Xl (s^) n or of order O 



cp < h ^ ho, the largest component of hgrc is 



xi(s + )n. 

3. if ^ h < h cp , the largest component of ^crc is either 
of order Inn or of order & X i(s + ) n - 

The item 3 extends to the values ^ h < h cp what hap- 
pens in item d) of the random cluster model for vanishing 
field. It implies that the thermodynamic and topological 
lines coincide there. This holds also at (fi cp ,h cp ), with a 
new behaviour in the sense that the giant component is of 
order n 2 / 3 or of order Q Xl r s+ )n. 

The Kcrtcsz line is given in this range of temperature 
(and field) by 



h { K } (fi) = h T (fi) for fi cp < fi £<«) or < h h cp 

(23) 

The item 2 shows that the giant component exhibits a 
jump, however, for given h > h cp , this giant component 
already appeared for lower values of fi. The topological 
transition line is determined, as in case q = 2, by eliminat- 
ing s between the mean field equation and the condition 



(17). This gives 



h%(fi) = 



9-1 



— In — — — for fi ^ ficp or h ^ h cp 



9-1 



(24) 

see Fig. 2. Below the Kertesz line h$ the largest compo- 
nent is of order Inn, while above hjp and outside of the 
segment (19), the largest component is of order 8 Xl ( So )7i. 
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Fig. 2: The thermodynamic transition curve: thick (solid, 
dashed) line; the topological transition curve: solid (thick, 
thin) line. 



The symmetry of the free energy is better seen if we 
parametrize the magnetization s by z G (±1) as follows: 

q-2 



q 



(25) 



2( 9 -l) 2( 9 -l) 

l+z 
2 ' 

, and the free energy decomposes 



The corresponding densities are given by x\ 
x 2 = ■■■ 



— X 9 ~ 2(q-l) 

immediately into its even and odd parts 



/(*) = 



-[(1 + z) ln(l + z) + (l-z) ln(l - z) 
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The first derivative of / yields the mean field equation (2) 
while the second derivative is 



/"(*) = 



1 



fiq 



1 



4(9-1) 



(27) 



Thus f is convex when fi ^ fi and has exactly one 
minima. When fi > fio it has at most two local minimas. 
The analyticity of so outside of /it is a consequence of 
f"(zo) > 0. The monotonicity of fix\(s+) is trivial, that 
of fixi(s-) can be infered from the computation of the 
contour lines for fix\(s) below /it- 
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Conclusion. — To summarize, we have given the 
equation of the Kertesz line of the mean-field Potts model 
and shown that when q 3 this line exhibits a cusp at 
some (/3 cp ,/i cp ). For low fields (h ^ h cp ) the Kertesz line 
coincides with the thermodynamic transition line, and for 
large fields (h > ho) only the topological transition re- 
mains while the thermodynamic transition disappears. In 
addition, at intermediate fields (h cp < h < ho) the Kertesz 
line separates from the thermodynamic line. This means 
that decreasing the temperature one sees first the appeare- 
ance of a giant component (on hx) and, then on hx, this 
component exhibits a jump. This behavior is new com- 
pared to what happens for the model on the 2-dimcnsional 
regular lattice. There, no intermediate regime appears [4]. 
We expect for sufficiently high lattice dimensions, a be- 
haviour similar to the model on the complete graph. 
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